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Abstract 

In this work we discuss the zero temperature limit of a "p-wave" holographic superconductor. The 
bulk description consists of a non-Abelian 5(7(2) gauge fields minimally coupled to gravity. We numer- 
ically construct the zero temperature solution which is the gravity dual of the superconducting ground 
state of the "p-wave" holographic superconductors. The solution is a smooth soliton with zero horizon 
size and shows an emergent conformal symmetry in the IR. We found the expected superconducting 
behavior. Using the near horizon analysis we show that the system has a "hard gap" for the relevant 
gauge field fluctuations. At zero temperature the real part of the conductivity is zero for an excitation 
frequency less than the gap frequency. This is in contrast with what has been observed in similar scalar- 
gravity-gauge systems (holographic superconductors). We also discuss the low but finite temperature 
behavior of our solution. 

1 Introduction 

In recent years string theory or, more specifically, gauge-gravity duality has seen interesting applications 
in the field of condensed matter physics. One of the earliest such applications is the discovery of a 
superconductor- like phase transition in AdS with a Reissner-Nordstrom black hole and a charged scalar 
field minimally coupled to a U(l) gauge field [DO- 

Earlier works on this system at low (but non zero) temperatures indicated the presence of a gap in 
the spectrum, with an exponentially decaying population of normal carriers following the Boltzmann 
distribution. In various recent works [3J [3] O [S] El El E] 5 authors have addressed the problem in the 
zero temperature limit. More precisely, they consider a charged scalar field and gauge fields in an AdSi 
background at zero temperature. The solution turns out to be a solitonic solution with a zero sized 
horizon. The authors find that the effective potential for small gauge field fluctuations vanishes near the 
black hole horizon. This implies that the normal component of the A.C. conductivity never vanishes, 
even at zero temperature, which in turn indicates that the superconductor is gapless. This is a new and 
interesting result different from what would have been naively expected. 

In this paper we consider non-Abelian gauge fields in the background of an extremal AdSi black hole. 
The setup is quite similar to so-called "p-wave" holographic superconductors [10] and known to have a 
superconducting transition^. We work at the zero temperature and found a fully gravity back-reacted 
solitonic solution to this system similar to [3J. Then we consider the frequency dependent conductivity 
by considering certain vector fluctuations on this background. This is an anisotropic system which shows 
different conductivity at different directions. We find the expected superconducting behavior for relevant 
currents. Interestingly we find that the corresponding effective potential does not vanish at the horizon (it 
goes to zero near the boundary) . Following the argument in [3j , we thus arrive at the conclusion that the 
holographic non-Abelian superconductor does have a finite gap for the relevant gauge field fluctuations. 

The plan of the paper is as follows, in section[2]we discuss the general setup and Einstein's equation for 
a non-Abelian gauge field coupled to gravity. In section [3] we discuss the instability of the RN black hole 
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solution. In section [3] we construct the zero temperature solitonic solution which is the superconducting 
ground state of the theory. We also discuss issues related to conductivity and show the existence of a 
hard gap. In section [5] we propose a possible form of the solution at non-zero temperature using the 
method of matched asymptotic expansion. 



2 Setup 

The Einstein- YM action for a non-Abelian gauge field with a negative cosmological constant is given by 

nn. 



C = j d^x^—g jFrF^j , (1) 

where _F M „ is the field strength of an SU(2) gauge field. 

To tally with our plan of considering anisotropy in the spatial direction, we choose the following 
ansatz for our metric 

ds 2 = -g(r)e- x(r) dt 2 + ^- + r 2 ( cirfdx 2 + dy 2 ) , (2) 
g(r) \ J 

and the gauge fields 0, 

A= A{r)r 3 dt + Bir^dx. (3) 
It is straightforward to find the non-zero components of the field strength 

F«„ = d»A a v - d„Al + iq[A„,A„] a (4) 

are 

Frx = —Fxr = B'(r), F? t = -Fl = A'(r), F 2 t = -F 2 x = qAB. (5) 
The energy momentum tensor is 

1 SSmatter 1 /1^2\ , 1 ot/3 rpa r^a 



and the non-zero components of the energy momentum tensor are 



o9^ \- A F + ^FM, (6) 



T rr = ~-^A' 2 + -^B' 2 + - f (qAB) 2 , 
4gc 2 4r 2 c 2 4<? 2 r 2 c 2 

T xx = \e-r 2 c 2 A' 2 + 9 -B' 2 -e-^^, 

T yy = \ e *r 2 A' 2 -JL B ' 2 + e- { ^f. (7) 
4 4c 2 4gc 2 



It is to be noted that here T xx ^ T yy due our anisotropic ansatfl 
The Maxwell's equations of A(r),B(r) are 



Af — ► A" + A' I * + \ + £- ) - ^f-A = 0, 



Ai^ B " + B'(^-^-^) + ^^B = 0. (8) 



6 Due to a repulsive term coming from the non-Abelian interactions, it is expected that a isotropic ansatz will have a quartic 
instability and would possibly have more free energy than the anisotropic ones 10 12] • We left the detailed discussion of these 
issues for a future study. 

7 HSH thanks Matthew Roberts and Gary Horowitz for discussions on this point. 
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The potential term for A(r) is different from the similar term in EYMH case due to the existence of a 
factor 4j. The diagonal Einstein equations give, 

1 c'\ /l 3c' c"\ „ e x ,, 2 g „, 2 „ q 2 A 2 B 2 
' 9 + — + - +3 = -A' 2 + -^-B' 2 + e^ 



r 2c J \r 2 rc c J 4 4r 2 c 4r 2 gc 

-t + c l(^ + 9'\ _ e*q 2 A 2 B 2 



g J g' 2 r 2 c 2 ' 

,(g' (2 X '\\ B' 2 x g 2 A 2 B 2 

The above equations are invariant under the following scaling symmetries: 

r — > air, (t, x, y) — > (t, x, y)/ai, g—>a 2 g, ^—+01^4, B — > aiB, (10) 
e x — » a\e x , t — » a2t, ^4^^/02. 
x^x/az, B — » a3_B. c — * 03c. 

The second scaling symmetry may be used to set x = at infinity and the third scaling symmetry may 
be used to set c = 1 at infinity, so that the asymptotic metric is that of AdSi. 
The fields have the following asymptotic behavior: 

A = »-£, B = B h + ^-, (11) 
r r 

where /i is the chemical potential and p is the charge density in the boundary theory. In what follows 
we will only consider the solutions for the field B which vanishes near the boundary, i.e. Bo — 0. 

3 Reissner-Nordstrom black hole solution and Instability 

To see when one expects hairy black holes at low temperature, one can study linearized perturbations of 
the extremal Reissner-Nordstrom AdS (RN-AdS) black hole 3 . The general RN-AdS solution is given 
by, 

c=l, X = B = 0, g = r 2 -±-(l + p 2 /4)+p 2 /4r 2 , A = p(l-l/r) (12) 
The temperature of the black hole (1121) is 



T-l&ZpFL^ ,13, 

For AdS-RN, T = (12 — p 2 )/16-7r. It is known that fluctuations of B x in this background develops a 
tachyonic mode [TD] at low temperaure and for sufficiently large q. The resulting solution is a super- 
conducting black hole with vector hairs. Following 0], we argue that such an instability persists at the 
extremal limit. At the extremal limit (T — 0), we have p = 2^3. The near-horizon limit of the extremal 
solution is AdS2 X R 2 with a metric, 

dr 2 

ds 2 = -6(r - l) 2 rit 2 + _ + dx 2 + dy 2 , A = 2V3(r-l) (14) 

Plugging this into the wave equation of B(r), and changing variables f = r — 1, we recover a wave 
equation for AdSi with a new effective mass, 

B,fF + \b,t — —z^-B = 0, m 2 // = -|-. (15) 

The instability to form SU(2) vector hair at low temperature then is just the instability of scalar fields be- 
low the Breitenlohner-Freedman (BF) bound for AdS2'- rn 2 BF — —1/4. Thus the condition for instability 
is 

q 2 > f. (16) 
It is reasonable to expect that such an instability takes the system to its superconducting ground state. 



3 




Figure 1: a vs q plot shows a almost linear relation with a slight dip for a for a value of q < 2.5. 



4 Zero temperature solution 



Being a single state without any degeneracy, a superconducting ground state does not have any entropy 
associated with it. This fact concur with the observations in 3,4,5 that the zero temperature limit of 
the superconducting black holes in EYMH theory have zero horizon size. Motivated by these facts we 
assume that same is true for a non-Abelian hairy black hole also and the ground state is a geometry with 
zero horizon size. We choose the following ansatz, near r — » + , (similar to that of EYMH case [3]) 



A~A (r), B ~ B - Bi(r), x ~ Xo - Xi 0)> + pi(r), 



c + ci(r). 



(17) 



All the terms with subscript 1 are sub-leading and go to zero, faster than the leading part where it is 
applicable, as r — > 0. Putting Eq. (|17[) in Eq. J5), we get find out the various terms in Eq. (|17|l . For 
example from the equation of motion for A, 



n 2 R 2 

I B 

A 



A 



„(-?) 



a = qBo/co, 



where we have used the observation Aq — > at the horizon. Following similar procedures, we get 



A ~ A e~ a/r , B~B Q (l 



e X0 q 2 A 2 _ 2a/ , 



4q 2 



CO 1 + 



X ~ Xo 



~2r~ 



' loa _2a/r 2 e xa A 2 a _ 2a/r 
-e , g ~ r ; e 



8r 2 
o 

4r 



(18) 



(19) 



where by rescaling (|10[) we set the coefficients Ao = 1 and xo = 0, Co = 1. After we solve the equations 
we again use the rescalings of g, c, \ to make the Asymptotic geometry the same as that of AdSi. For a 
given q, we choose a in such a fashion that B vanishes near the boundary jjj . The numerical evaluation 
shows us an almost linear relation between q and a (Fig[T]),with a slight dip for a lower value of q. Fitting 
a linear relation between q and a with a range of 3 < q < 8, we get a w 0.088 + 0.66g. 

Our numerics is less reliable in the region q < 0.95. It is not clear what exactly happens in this 
region. As q — > w 0.866, the geometry of our solution seems to approaches extremal RN geometry 
(Fig[T]). However the questions of phase transitions remain unclear. It is expected that there is a second 
order phase transition at q — ^ from RN black hole to our solution. However one may not rule out 
the possibility of a first or any other type of phase transition(s). At least in the range of our numerics it 
seems that our solutions have lower free energy than the extremal RN black hole in a grand canonical 
ensemble. As expected the solitonic solution is more massive than the RN black hole, but they carry 
more charge for the same boundary chemical potential and consequently have a lower free energy. 

General properties of our solution is similar to EYMH case [3]- B,x approaches a constant near 
r = 0, g — r 2 and the gauge field A and its derivative (electric field) vanishes. The metric approaches 
AdSi with the same value of the cosmological constant at infinity. The extremal horizon is just the 
Poincare horizon of near horizon geometry AdSi. In terms of the dual field theory, this means that the 
full conformal symmetry is restored in the infrared. Due to the asymmetry of our solution in x and y 
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Figure 2: The left panel shows a plot of the rescaled with the rescaled variable r n . From top to bottom 
q = 1.5, 1.1, 1, 0.9. The dashed line is for the RN black hole. The right panel shows a plot of A with r n for 
the same set of values of q, from bottom to top. 



direction the ratio of speed of light in IR and UV is not only non-unity but also different in different 
directions. We have, 

^ = ^ exp( T ) ' <^ = exp( T ) (20) 

Our solution may be viewed as static, charged solitons, where the electrostatic repulsion is balancing 
the gravitational and non-Abelian interactions. As can be checked easily from the asymptotic behaviors 
the energy of the soliton is finite. One interesting feature is the existence of an "essential singularity" (i.e. 
an e~~ behavior) of A(r) at r = 0. Similar singularities appear for other quantities in Eq. (|20[1 . Hence 
all the sufficiently high derivative of the metric and gauge helds vanishes at r = 0. This is in contrast 
with the EYMH case [3| , where sufficiently large curvature invariants diverge at r — for a generic a and 
behavior of the helds near r = is power law. In terms of the "tortoise" co-ordinates (|27[1 , the essential 
singularity is just the reflection of the fact that A is massive near the horizon and decays exponentially. 

4.1 Conductivity 

In order to calculate the conductivity of this system, we need to turn on a small perturbation in the 
vector potential. The gauge field and the associated metric perturbations are of the form: 

Ay = ea(r)e _u "V 3 di/, g ty = eh(r)e~^\ (21) 

where we assume that both the gauge and the metric perturbations are of the same order 0(e). We won't 
be looking in the other component of the gauge field perturbations which decouples from the above. From 
Maxwell's equation we get, 



y V Q 2 c / V Q 2 gr 2 c 2 ) g \ \ h 2 c / r 2 gc 2 



(22) 



To linear order, Einstein's equations give, 



-h + h' = -aA', (23) 



r 



8 However it should be noted that B(r) is not square integrable. This is not a requirement of a valid solution to the Einstein's 
equations and assuming that our numerics are reliable, such solutions exists. The question is that whether they play a role in 
the context of AdS/CFT. In a sense whether it belongs to same quantum ensemble as the extremal black hole. As our zero 
temperature solution is of finite energy this seems to be the case. There are well known solutions such as the BPST instanton 
|19j which also have a non square integrable gauge field configuration with a bounded action. In our case, just as in case of 
BPST instanton, the solution may be approached by configurations which are of finite norm with a bounded energy (just replace 
B — » cxp(— ^)B(r), /3 — > 0). At finite temperature the gauge field B is likely to be square integrable due to the cutoff provided 
by the horizon. Our zero temperature solution may also be thought as a zero temperature limiting case of such a solution (see 
Sec. 0. 
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Figure 3: Plot of the rescaled potential V(r n ) for q = 1.5, 1.1, 1 (from the right) at T = with rescaled 
co-ordinate r n . The poincare horizon is at r n — for the zero temperature solutions. V(r n ) is non-zero at 
the horizon. The dashed curve is for the extremal black hole with a horizon at r = 1. 



which combined with the equation of A(r) to give: 

*" + *'(--4 + -)+ a ( e -^- :L ^- eX — ) =°- ,2J! 
\g 2 cj \ g 2 gr 2 c 2 g ) 

This can be written as a Schrodinger equation: 

— a" + V(f)a = <?uj 2 a, (25) 

where 



V(r)=g(SA' 2 + e-*^). (26) 



and all the derivatives in Eq. (|25(l are in terms of the new variable new variable f ( "tortoise coordinate" ) 
given by: 

^ — -x/2 « fr,r,\ 

d~r= e 9C Tr- (27) 
Here, c approaches to unity as r — » oo, so that the spacetime is asymptotically yldS'4. It follows then 
from Eq. that the potential V(r) vanishes near the boundary. If we require g ~ r 2 near the horizon 
at r = then the first term vanishes, while the second term is finite as B(r = 0) = Bo 7^ and \ is also 
finite at the horizon. The value of the potential at the horizon decreases as q — > ^ (Fig(3|. Note that 
since c — > 1 near the boundary, the quantity ui can be interpreted as the frequency of the incoming wave. 

The superconducting nature of the system is argued from the existence of a supercurrent solution. If 
we choose w = and integrate A% from the horizon (with a regularity condition at the horizon), we are 
expected to get a non-trivial A x . Existence of such a solution implies a 8 function for the real part of 
conductivity at w = [20ll3]Fl 

Re(<r(w)) ~ 5{uj) + finite. (28) 

The nature of the finite part of the conductivity can be inferred from the potential V(r). As argued 
before the potential vanishes near the boundary and is finite near the horizon. The fact that the potential 
is nonzero at the horizon makes it possible for this system to exhibit a hard gap. From Eq. (|25|l . the 
field a has the following asymptotic behaviours near the horizon (f — » —00) and the boundary (f — ► 0): 

a(f -> 0) ~ a + a\f (29) 

(30) 



9 One should be more careful about the existence of a delta function for the real part of conductivity at cj = and relation 
between superconductivity. Due to the coupling of gravity modes the potential V(r) is non-trivial even in the normal phase 
(Fig [3} and consequently a delta function is present in the real part of DC conductivity. Physical interpretation of such a delta 
function comes from the "translational invariance" of the theory. See 12 1 1 for further discussions. 



G 



where Cj = \J CqLu 2 — Vq, with c; t , being the near-horizon values of c and a respectively. Here, we have 
chosen the incoming boundary condition near the horizon. The conductivity of the system is given by: 



id 



a — — 



(31) 



It follows from Eq. (|25|) that : 

a* a" — aa* — 0, (32) 

which implies that the quantity A = a*a! — aa* = 2ilm(aa* ) is a constant. Equating the values of A 
near the horizon and the boundary we get: 

£ \ a o\ 2 ,~,2 v. n 

Re(er) = <{ -k^F ' U >U (33) 

0, <Ir < 



Therefore, the real part of the conductivity will vanish whenever Cj is imaginary, i.e. when u < yVo/co, 
which defines the gap. 

As our solution is an anisotropic superconductor (only superconducting in A y direction), the ex- 
citations in A^. channel may not be gapped [TO]. The gap seems to be a property of the anisotropic 
non-Abelian system, as the scalar-Abelian gauge field system does not exhibit this gapped behaviour in 
the extremal limit [3]- 



5 Finite temperature analysis 

Although the Eq. (JHJ has been solved in a probe limit, the full set of equations with gravity backreaction 
(Eq. ©) is not yet solved at a non-zero temperature. Hence it is an issue whether our zero temperature 
solution may be realized as a zero temperature limit of the finite temperature hairy black hole (and 
whether such a low temperature hairy black hole solution at all exists.). One affirmative clue comes 
from the emergent conformal symmetry in the IR. In the near horizon region, i.e. a 3> r, our solution 
approaches a AdSi geometry with a background gauge field A x = Bo and At = 0. One may consider 
a black hole situated deep inside this AdS space, such that the radius of the black hole ro <C a. Here, 
g(r) near the horizon is given by, g(r) = r 2 (l — However this Schwarzschild like solution will get 
correction due to a small non-zero A(r) near the horizon. The solution of A(r) in the r <C a is given by, 

A(r)«exp(-£)(1-^). (34) 

A(r) is small in the region r<a. In this region we can do a perturbation in A(r) for other quantities. 
Realizing that our scaling relations Eq. (|20|) is just a perturbation in A(r), we can match a near horizon 
solution with the scaling solution (i.e. Eq. (|20[l ) in a region ro C r < a. Hence we use the method of 
matched asymptotic expansion to argue that, 





A T (r) = 


A(r)(l- 


ro 
r 


B T {r) = 


Bo - (B(r) - 


Bo)(l- 


Toy 

r 


c T (r) 


= c + (c(r) - 


-co)(l- 


TO: 

r 


Xr{r) = 


X(r) + (x(r) 


-Xo)(l- 


ro 
r 




ffr(r) = 


= .9(0(1- 


ro 
r 



(35) 



is a solution to Einstein's equations if ro <S a. The quantities (A, B.c, x) appear in the right hand 
side of the above equations are our numerical solution and the quantities with a subscript T are the 
finite temperature versions. Whether correction to such a solution is finite or not is an open question, 
which may be partially settled by a convincing finite temperature numerics. However due to the possible 
"smooth" nature of the equations Eq. @, i.e. a small change in the near horizon boundary condition 
leads to a small change near the boundary, it is reasonable to believe that corrections to Eq. (|36p are 
well controlled. What happens to the potential V(r) (see Eq. (|26[) ) in this limit is interesting. At any 
non-zero ro, the potential V(r) goes to zero at the horizon (r = ro) and the height of the potential 
is always finite. Hence the transmission through this potential is non-zero. In our language the (i.e. 
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Figure 4: Schematic plot of the potential V(f) for different T. T decreases from the left. V(r) is gradually 
widened as T decreased. Horizon is at f = oo. 



Eq. (J33J ) u — uj is a real number. Hence the real part of the conductivity is non-zero for any value of u). 
However, in the "tortoise co-ordinates" (|27p the potential is gradually widened as vq — > (FigfiJ). Hence 
the transmission through the potential gradually decreases as the horizon size is decreased. We expect a 
exp(— fall off for the real part of conductivity. A detailed discussion of the various gaps etc is left for 
a future work. 
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